Bloch points are sticky 
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Bloch points are magnetic topological defects resembling hedgehogs. The discrete nature of a 
magnetic lattice creates a periodic potential that can pin a Bloch point. The pinning force is of the 
order of the exchange constant, a few piconewtons in a typical ferromagnet (permalloy). A domain 
wall containing a Bloch point can have a sizable depinning field in the tens of oersted. 

I. INTRODUCTION 

Topological defects, nonlinear excitations made stable by their nontrivial topology, are ubiquitous in magnetism. 
They play an important role in the dynamics of magnets, mediating transitions between low-energy states. Defects 
come in a variety of dimensions and shapes. In a thin ferromagnetic film, where shape anisotropy keeps magnetization 
in the plane of the film, point-like topological defects are vortices. Away from the core, magnetization stays in the 
easy plane, winding clockwise or counterclockwise once as the core is encircled. At the core of a vortex, an area several 
nanometers across, magnetization points normal to the easy plane. The core endows the magnetization with a double 
twist, which is ultimately responsible for a gyrotropic force that dominates vortex dynamics. 

An example of a point-like defect in three dimensions is a magnetic hedgehog, or a Bloch pointer—. In the simplest 
case, magnetization in its vicinity has a spherically symmetric distribution, 
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where M is magnetization length, with a typical value of 8 x 10 5 A/m (permalloy). More generally, a hedgehog is 
stable if the unit vector of magnetization rii(r) = M(r)/M has a nonzero winding number of the second homotopy 
group TT2(S 2 )— also known as the skyrmion charge: 
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where the integral is over a closed surface surrounding the hedgehog. For a simple hedgehog dU, q = ±1. 

A skyrmion number q can also be defined for a vortex on a two-dimensional surface. The skyrmion density 
eijk m ■ djrii x (9 fc m/87r is nonzero only if the unit vector of magnetization m is non-coplanar. Therefore, the skyrmion 
charge of a vortex is concentrated near its core, where magnetization comes out of the easy plane. The skyrmion 
charge of a vortex is half-integer, q = ±1/2. This difference between a hedgehog and a vortex can be understood 
intuitively by noting that magnetization m around a hedgehog acquires every direction on the unit sphere, whereas 
magnetization near a vortex core only covers half of it, e.g., the northern hemisphere. The skyrmion number of a 
vortex on the surface of a magnetic sample is related to its "core polarization," i.e., the direction of magnetization 
normal to the surface. For a vortex, a core magnetized into (out of) the surface has q = +1/2 (—1/2). For an 
antivortex, the relation is reversed. 

Whereas magnetic vortices have been extensively studied both theoretically and experimentally and are by now well 
understood, isolated magnetic hedgehogs are elusive and have so far been studied mostly by theorists. (Kabanov et ai— 
observed Bloch points on a Bloch line, which in turn is part of a Bloch domain wall.) These objects are fascinating 
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FIG. 1. (a) The transverse wall in a cylindrical nanowire of radius 10 (nm). The left (right) figure is a snapshop of magnetization 
at the center of the cylinder from the positive x (z) direction. The blue (red) arrows mean the positive (negative) component of 
magnetization along the x axis, (b) The Bloch point wall in a cylindrical nanowire of radius 30 (nm). The same color scheme 
as (a) is applied. 
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as well as important. For instance, the reversal of core magnetization in a vortex changes its skyrmion number by 
±1. One might suspect that the difference is taken by an object with unit skyrmion charge such as a Bloch point. 
Indeed, numerical studies have demonstrated the involvement of a Bloch point in a core reversal process.— ^ Hedgehogs 
also create nontrivial magnetization dynamics. Hertel and Kirschner— reported an observation of "magnetic drops" 
behind a moving Bloch point. 

Another peculiar feature of a magnetic hedgehog is the absence of a characteristic length scale. Whereas the size 
of a vortex core is determined by a competition between short-exchange force and shape anisotropy due to dipolar 
interactions, a hedgehog solution (fTJ) is scale-invariant and thus sensitive to both short and long-wavelength physics. 
We may expect that the center of a hedgehog, where spatial variation of magnetization is particularly strong, can be 
easily pinned by lattice imperfections and even by the atomic lattice itself! In this paper we show that this is indeed 
the case: even the periodic potential of an atomic lattice may create a substantial landscape for a Bloch point. 

We present a theoretical study of a Bloch point that is part of a domain wall in a ferromagnetic nanowire with a 
solid circular cross section (a solid cylinder). Thanks to shape anisotropy of dipolar interactions, such a wire has two 
ferromagnetic ground states with uniform magnetization along the axis of the wire. A domain wall separating two 
ground states can be of different types, depending on the radius R of the wire^ - — . For a small radius, not exceeding 
a few exchange lengths, a "transverse" wall is formed, Fig. HJa). It consists of two surface defects, a vortex and an 
antivortex with skyrmion numbers q — ±1/2 each (negative for a head-to-head domain wall). For a larger radius, a 
"vortex" wall is formed, Fig.[T](b). Contrary to its name, it contains a Bloch point, rather than a vortex, so we shall 
refer to it as a Bloch-point wall (BPW). The skyrmion charge of the hedgehog is ±1 (negative for a head-to- head 
domain wall). 

The domain-wall topology of magnetization in a magnetic wire provides an easy way to apply a force to a Bloch 
point: a uniform magnetic field Ho parallel to the axis of the wire pushes the domain wall with a force [IqQHq, where 
Q = ±2itR 2 A4 is the magnetic charge of the domain wall (positive for a head-to-head domain wall). If the Bloch 
point is pinned, the wall is unable to move until the applied field is sufficiently strong to overcome the pinning force. 

The pinning of a Bloch point in an atomic lattice has a simple origin. Picture a hedgehog moving through a simple 
cubic lattice of magnetic atoms parallel to one of the cubic axes, Fig. [5] The hedgehog's energy is lowest when it 
is at the center of a cubic cell (z — 0), farthest from the magnetic dipoles in its corners. It is highest (along that 
trajectory) when the hedgehog is at the center of a face (z — a/2). A simple model of the Bloch point's energy on 
this straight line is U{z) — — Uq cos (2ivz/a), where a is the lattice constant. The amplitude Uq can be estimated as 
follows. The primary source of the pinning potential is the exchange energy E CK — A j cPrdim ■ dim, where A is 
the exchange constant with the dimension of energy per unit length. We thus expect that Uq = cAa, where c is a 
numerical constant of order 1. This gives the maximal pinning force F p = 2ttcA that is independent of the lattice 
constant a. The critical field for depinning a domain wall with a Bloch point in a magnetic wire of radius R is then 
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where £ cx = \j2Aj (j,qM 2 is the exchange length. For a typical ferromagnet (permalloy, A = 1.3 x 10~ n J/m, 
M = 8 x 10 5 A/m, £ cx — 5.7 nm), we obtain a pinning force of order 10~ n N. A permalloy wire with radius R = 50 
nm would have a sizable depinning field H c « 60 Oe. 

The paper is organized as follows. In Sec. [Til we compute the lattice potential of a magnetic hedgehog and the 
associated critical field for a Bloch-point domain wall. In Sec. IIII1 we compare the computed value against numerical 
simulations using the micromagnetic simulator OOMMF.— We discuss our results in Sec. IIVI 



II. THEORY FOR THE PINNING CRITICAL FIELD 



The energy of a ferromagnet in a continuum model is given by 

E = J d?r (a d t m ■ d. t m - ^-m ■ h D - Km"*) , (4) 

where m = M/M is the unit magnetization vector, A is the exchange constant, M is the saturation magnetization, 
K is an easy-axis anisotropy constant, and ho is a dimensionless demagnetizing field satisfying V • ho = —V ■ m 
and V x hp = 0. The coupling constants A, /ioM 2 , and K have dimensions of energy/length, energy/length 3 , and 
energy/length 3 , respectively. They can be used to construct an intrinsic length scale known as the exchange length, 
£ cx = \j2Aj ^qM 2 , and a dimensionless "quality factor," Q = 2K/fi M 2 expressing the strength of local anisotropy 
relative to the shape anisotropy of the demagnetizing field hjj . Since we are interested in a weakly anisotropic system, 
Q <C 1, we ignore the anisotropy energy for calculation of the mesh potential of BP. 
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FIG. 2. Magnetization configurations with two different BP positions (cross): (a) at the center of a cubic cell and (b) at the 
center of a cell's face. 



We have evaluated the energy of a magnetic hedgehog in a Heisenberg ferromagnet on a simple cubic lattice with 
nearest-neighbor exchange coupling J — 2Aa, which gives the same exchange energy for slow spatial variations as 
does the continuum version, Eq. The lattice Hamiltonian is 



E = -2AaJ2 
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where is the unit vector pointing from the site i to the site j and r.y is the distance between two sites in units of 
the mesh size a. Spins are located at half-integer lattice points = (ii + 1/2, «2 + 1/2, «3 + 1/2), where i\, 12, and iz 
are integers. We used the following Ansatz for magnetization in cylindrical coordinates r = (jj, <p,z), z being the wire 
axis: 
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where (0, 0, zbp) is the location of BP, (f>o is the deviation magnetization projected into the xy plane from the azimuthal 
direction (p, and s is a phenomenological shape factor expressing the deformation of a hedgehog inside a cylindrical 
wire. The shape factor was determined from numerical simulations by fitting magnetization of a static hedgehog to 
Ansatz in the vicinity of the singularity. 

Our analysis, presented in detail in Appendix shows that the periodic potential experienced by the hedgehog in 
the lattice is dominated by the first harmonic: 



E(zbp) « -Aa(c cx - c dip a 2 /€g X ) cos (27rz B p/a), 



(7) 



where c cx and Cdi P are positive constants of order of 1 reflecting the contributions of exchange and dipolar interactions. 
Figure [3] shows the dependence of c ex and Cdi p on the shape factor s. In the case of the atomic lattice, a <C i cx and the 
dipolar term contributes relatively little. However, in our numerical simulations the lattice constant a was comparable 
to the exchange length £ cx . We therefore kept the dipolar term in Eq. j7} for comparison with the numerics. 
The pinning force created by the periodic potential is 



Fp(z BP ) = -E'(z BF ) w 2nA(c, 
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It matches the external force of the applied field, Fz = poQH, where Q = 2itMR 2 is the wall's magnetic charge, until 
F p reaches its maximum value at zbp = a/ A. This gives a critical field for depinning, 
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III. NUMERICAL SIMULATION 



We performed numerical simulation with the micromagnetic package OOMMF 16 to verify our calculation of critical 
fields. The simulated sample has the following material parameters: exchange constant A = 1.0 x 10~ n J/m, 



4 



0.5* 

0.4 
0.3 
0.2 
0.1 



■ c dip 



0.6 0.8 1.0 1.2 1.4 

FIG. 3. Dependence of c cx and Cdi p on the shape factor s. 
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FIG. 4. (a) The position of the Bloch point and external field as a function of time (b) The velocity of the Bloch point as a 
function of its position 

saturation magnetization M — 10 6 (A/m), easy-axis anisotropy K — 10 4 J/m 3 , gyromagnetic ratio 70 = 2.21 x 10 5 
m/(A s), and a damping constant a — 0.3. (We used a large damping constant to reduce simulation time.) The 
exchange length for this material is £ cx — 4.0 nm and the quality factor is Q = 0.016. A small easy axis anisotropy 
helps stabilize a domain wall with a Bloch point. 1? Since Q <C 1, we did not include the anisotropy energy in the 
calculation of the mesh potential. We used two different sample geometries: (i) radius R = 10 nm and length L = 200 
nm; (ii) radius R — 15 nm and length L — 300 nm. To determine the critical field, we applied the external field 
increasingly starting at 100 Oe and 30 Oe and increasing it in steps of 10 Oe and 5 Oe for the narrow and wide wires, 
respectively. Figure 2] shows the position of the Bloch point zbp as a function of time for the narrow wire with a 
lattice constant a = 1.7 nm. 

The results are summarized in Fig. [5] It can be seen that the depinning field H c does not extrapolate to zero as 
the size of a unit cell a in the simulation is reduced. This behavior of a Bloch-point domain wall is strikingly different 
from the typical behavior of domain walls in numerical simulations. Ordinarily, the effects of lattice discretization 
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FIG. 5. Theoretical and simulated critical fields 
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become unnoticeable when the mesh size is reduced below half the exchange length £ cx . The hnest mesh size in our 
simulations, a = 1 nm, is much smaller than £ cx = 4.0 nm, yet the lattice effects do not subside. Furthermore, the 
critical held increases as the mesh size decreases. The main reason for this effect comes from the dependence of the 
shape factor on the mesh size. As the mesh size decreases, the shape factor decreases, which results in increase in c ex 
(see the figure [3]). Magnetization around BP becomes more isotropic as the mesh size decreases. 

The mesh potential also affects dynamics of a BPW under the external field above the critical field. The velocity of 
BPW oscillates and this oscillation is periodic in the position of BPW, which means that the mesh potential for BP 
is responsible for the oscillation. See (b) in the figure 0] for an example. The slight increase of the average velocity 
as a function of time can be attributed to the attractive force between the volume magnetic charge of BPW and the 
surface magnetic charge at the end of the cylinder. 



IV. DISCUSSION 



Topological defects in magnets can be pinned by a periodic potential of the underlying magnetic lattice. When 
lattice discretization is employed to simulate a continuous magnetic medium, this periodic potential is an artifact that 
can be removed by choosing a sufficiently fine mesh. For a domain wall, the lattice period should be small compared 
to the wall width; for a vortex, smaller than its core. A magnetic hedgehog turns out to be exceptional in this regard 
because it has no characteristic length scale. As a result of that, the amplitude of the periodic pinning potential of the 
lattice is linearly proportional to the lattice constant. Thus the maximal pinning force is independent of the lattice 
constant and is of the order of the exchange constant, A = 1 .3 x fCP 11 N in permalloy. 

In this paper we have estimated the strength of the periodic pinning potential experienced by a magnetic hedgehog in 
a simple cubic lattice. The estimate is in good agreement with micromagnetic simulations, which yield a finite pinning 
field for a domain wall containing a Bloch point, even when the unit cell is much smaller than the exchange length. We 
expect that Bloch points are strongly pinned by the atomic lattice in ferromagnets. The strength of pinning should 
be reduced by thermal and quantum fluctuations, which can create a core where magnetization vanishes. Numerical 
cstimates&il show that the core size induced by thermal fluctuations is small, comparable to the lattice spacing at 
room temperature. 
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Appendix A: Mesh potential calculation 

The exchange and dipolar energies of a ferromagnet (O in the presence of a Bloch point at (0, 0, zbp) are 

E cx (z BP ) = -2Aa^rhi • m 3 

m 

EaUzbp) = kT^TF ' (A1) 

with rhi given by Ansatz ([5]). Being periodic functions of zbp with a lattice period a, these expressions can be written 
as Fourier series: 



E cx (z B p) - E cx (0) = -Aa ^ c cx,« cos(27rnz B p), 
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The subtraction of the zbp = term improves convergence of the energy in a finite lattice. The Fourier amplitudes 
c GX ,n and Cdip, n depend on the shape factor s; Cdi P , n also depends on the angle 4>q. Overall signs of formula is determined 
so that the first harmonics c cx l and Cdi P; i are positive. The shape factor s and the angle 4>q were obtained by fitting 
magnetization of an OOMMF solution at the eight mesh points closest to a Bloch point at zbp = 0. 
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The Fourier amplitudes c ex , n can be computed analytically by applying Poisson's formula to the index 13 along the 
z-axis, 
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We thus obtain 
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Here jq = (ii + 1/2, + 1/2,0) is the projection of r; onto the xy plane and <fo = tt/2 is used. This expression is 
dominated by the n = 1 term. In practice it suffices to use only n = 1 and «i, Z2 = — 1, 0, 1. c ox used in the main part 
of the paper is c ex ,i calculated with ii,Z2 = — 1,0,1. 

The Fourier amplitudes c cx .i and Cdi Pl i have also been evaluated numerically. After setting the maximum value on 
/ = max(|ii|, |i2|, in Eq.f IA1[) . we computed the dependence of E cyL and -Edip on zbp by increasing zbp from to 
1 in steps of 0.1. We then extracted c 0Xi i and Cdi P ,i by a Fourier transformation. These coefficients converge rapidly 
as I increases and I = 10 provides good accuracy. Cdi p used in the main part of the paper is numerically calculated 
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